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Vortex structures in soft magnetic nanodisks are highly attractive due to their scientific beauty and potential 
technological applications. Here, we experimentally demonstrated the resonant amplification of vortex 
oscillations by application of simple coherent field pulses tuned to optimal width and time intervals. In order 
to investigate vortex excitations on the sub-ns time scale, we employed state-of-the-art time-resolved 
full-field soft X-ray microscopy of 70 ps temporal and 25 nm lateral resolution. We found that, due to the 
resonant enhancement of the vortex gyration motion, the signal input power can be significantly reduced to 
~ 1 Oe in field strength, while increasing signal gains, by increasing the number of the optimal field pulses. 
We identified the origin of this behavior as the forced resonant amplification of vortex gyration. This work 
represents an important milestone towards the potential implementation of vortex oscillations in future 
magnetic vortex devices. 



P 



ast achievements in spintronics have largely relied on an endowment of electron-charge-based electronics 
with novel magnetization excitation phenomena such as single-domain switching 1 ' 2 , vortex-core switching 
in nanomagnets 3 " 18 , domain-wall motions in nanowires 19 " 23 , and spin-wave excitations in magnonic crys- 
tals 24,25 . Ultrafast (~ps) magnetization excitations and their collective modes in nanomagnets have attracted 
interest owing to their suitability as information carriers with the additional benefit of much reduced energy 
dissipation compared to other electronics-based circuit devices. 

A prototypic example is the study of the fundamental dynamics of vortex gyrations in single nano-size disks 
and their coupling in one- or two-dimensional arrays of vortex-state disks 26 " 29 . Since the gyration mode, repre- 
senting the orbiting motions of an individual vortex core around its equilibrium position at characteristic 
frequencies from several hundred MHz up to ~1 GHz 30 " 34 , is well understood and lends itself to accurate 
analytical predictions 35 , the vortex gyration has attracted much attention it can be used as carrier for informa- 
tion 26-29 or as a source for microwave emission 36 ' 37 . For implementation in future potential devices, however, the 
gyration amplitude must be large enough for sufficient signal gains and signal transfer fidelity in integrated 
functional blocks along with significant reductions in the operating power consumption 26 ' 27,38 . 

Here, we experimentally investigate the resonant amplification of vortex oscillations in nanodisks by optimally 
tuned field pulses, by time-resolved soft X-ray microscopy measurements. The application of a series of coherent 
magnetic field pulses of ~ 1 Oe triggers a resonant enhancement of the gyration and subsequently an increased 
vortex-core oscillation amplitude and, therefore, signal gain. This mechanism provides a novel and efficient way 
to operate reliably future devices at significantly reduced power consumption. We show that this behavior can be 
explained by individual force terms and their balance, as obtained from combined results of micromagnetic 
numerical and analytical calculations. 

Results 

Vortex oscillation. Figure la shows the spin configuration, i.e., the magnetic vortex structure, in a single circular 
magnetic disk. It is characterized by an out-of-plane magnetization (referred to as polarization) in the core region 
which can point either upward or downward, and an in-plane curling magnetization (referred to as chirality) 
around core which can point either counter-clockwise or clockwise 39 " 41 . The vortex is a stable ground state and 
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Figure 1 | Concept of resonant amplification of vortex oscillations by 
coherent serial pulses, (a) Top: Perspective -view image and cross-section 
of magnetic vortex state. The color and height of the surface indicate the in- 
plane and out-of-plane components of local magnetizations, respectively. 
Bottom: Potential energy surface and analytically calcualted vortex-core 
trajectories. The vortex-core trajectories under the Gaussian-pulse train 
(pulse number N = 2) are plotted on the potential energy surface. For the 
cases I and II, the time intervals employed were AT = 2.0 x (2iz/cdd) and 
AT = 1.5 x (2k/cqd), respectively. The blue and red lines indicate the 
trajectories of vortex-core positions driven by the first and second pulses, 
respectively (see the pulse profiles in (b)). (b) Profiles of train of Gaussian- 
pulse currents along single-strip electrode with pulse width (er), time 
interval (AT), and pulse number (N = 4). 

results from the competition between the short-range exchange and 
long-range magnetostatic interaction. For small deviations of a 
vortex core from its equilibrium position, the motion of a vortex 
core with its velocity v is known to be well described by Thiele's 
equation 42 ' 43 : -G x \-Dv + dW(X)/dX = 0, where G=-pGz is 
the gyroforce vector of a given core polarization p, and D = DI is 
the damping tensor with the identity matrix I and its damping 
constant D < 0 [Ref. 35]. The W(X) is the potential energy for a 
displaced core position X and can be expressed as W(X) = 

W(0) + -k;|X| 2 + Wh> where k is the stiffness coefficient, the first 

term W(0) the potential energy for a vortex at its initial position X 
= 0, and the second term the potential energy for the shifted vortex- 
core position 35 . The last term, W H , is the Zeeman energy due to a 
driving force if any field is externally applied. With this approach, the 
force-balance equation is essential to understand and predict vortex 
motions, which is given as F G + F D + F w + F H = 0, where each term 
corresponds to the gyroforce F G , the damping force F D , the restoring 
force F w , and the Zeeman force F H , respectively 35,42 ' 43 . If there is no 
external driving force and damping, the restoring force F w is in 
balance with the gyroforce F G , resulting in a persistent motion of 
the vortex core. If, however, there is intrinsic damping of the 
constituent magnetic materials (nonzero damping parameter a) a 
damping of core oscillations will result. In order to allow vortex- 
core motions to compensate or overcome such damping, external 
forces such as the static, alternating, or pulse form of magnetic fields 
or spin-polarized currents need to be applied. 

Tuning resonant amplification of vortex gyrations. Here, we focus 
on a mathematically simple but technologically very useful pulse 
form as a driving force. Figure lb shows a serial -pulse train 
composed of Gaussian-shaped unipolar pulses of identical height 
I 0 and width o and time interval AT. In the case of a damped 
oscillator, the motion of the vortex core can be resonantly 
amplified by applying a series of driving forces synchronized with 



the oscillatory motion of the vortex. Arbitrary driving forces would 
increase the oscillation amplitude whenever the driving force is in the 
same direction as the motion of the oscillator, but would result in a 
forced damping when in the opposite direction. Thus, to resonantly 
build up vortex-core oscillations, the driving force needs to form 
an in-phase repeating pattern of the same frequency as the 
eigenfrequency of the oscillator. Similarly, if a train of Gaussian 
pulses of width (a) and time interval (AT) (as shown in Fig. lb), 
tuned to the intrinsic resonant frequency (w D ) of the given vortex, 
were applied either in-phase or out-of-phase relative to the vortex 
gyration, the serial successive forces would result in acceleration or 
damping of the vortex gyration, according to the time interval. The 
trajectories of the vortex core driven by two Gaussian pulses are 
illustrated in Fig. la. When the time interval is tuned to an integer 
of the vortex' precessional period, the amplitude of the vortex 
gyration can be resonantly increased, as seen in case I. By contrast, 
the vortex gyrations are rapidly decreased if the time interval between 
the first and second pulses is a half integer of the vortex precessional 
period, as seen in case II. On the basis of this approach, we used 
unipolar field pulses with an optimal width a = l/co D and coherent 
time interval AT = 2tz/co d (see Supplementary text online). 

Soft X-ray imaging of vortex oscillations. In order to 
experimentally demonstrate the above described concept of 
resonant vortex oscillations, we fabricated a sample, as shown in 
Fig. 2. The sample contains several single dots, each of which 
consists of Permalloy (Py: Ni 80 Fe 2 o) of diameter 2R ranging from 
2.0 to 6.0 /mi and a given thickness L = 40 nm (for experimental 
details see Methods). Trains of Gaussian current pulses applied along 

N-l 

the y axis are simply expressed mathematically by l y = I 0 ^ exp 

( — (t — mAT — At) /(2er 2 )), where N is the number of pulses in a 
given train. Core motions in the Py disk of diameter 2R = 3.0 fim 
were obtained by spatiotemporal- resolved full-field magnetic 
transmission soft X-ray microscopy (MTXM) which features a 
70 ps temporal resolution and a 25 nm lateral resolution utilizing 
a stroboscopic pump-and-probe microscopy technique 44 . X-ray 
magnetic circular dichroism (XMCD) contrast near the Fe L 3 
absorption edge (around 707 eV) provides magnetization orien- 
tation contrast. In our experimental setup, the in-plane curling 
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Figure 2 | Sample geometry for time-resolved soft X-ray microscopy 
measurement, (a) Schematic drawing of sample geometry consisting of Py 
disks with indicated diameter 2R = 3.0 /urn and thickness L = 40 nm and 
electrode of Ti (5 nm)/Cu (70 nm)/Au (5 nm). The electrode is positioned 
on the /-axis, thus yielding magnetic fields along the x-axis. The left-upper 
inset shows an optical microscopy image of the sample wherein Py disks of 
different diameters 2R = 2.0—6.0 /mi are placed. The sample surface was 
mounted at an angle of 60° with respect to the direction of the incident X- 
rays. (b) MTXM image of ground vortex state of Py disk of diameter 2R = 
3.0 fim and thickness L = 40 nm, obtained at the Fe L 3 edge. The vortex 
state is characterized by the CCW in-plane and downward core 
orientations. 
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magnetizations of the vortex were directly imaged as seen in Fig. 2b. 
The polarization and chirality of the Py dot were downward (p = — 1) 
and counter-clockwise (C = +1), respectively, as determined by the 
rotation sense of the vortex gyration and the in-plane MTXM 
contrast. The optimal pulse- width (er) value, as obtained from an 
experiment with single Gaussian -pulse excitation (N = 1), was 
1.45 ns. From this value, we could estimate a possible 
eigenfrequency of the real sample, here w D = 2n X 110 MHz, 
according to the explicit analytical form of a = l/w D [Refs. 17, 18]. 

As an illustration of the nontrivial motion of the vortex core in the 
experiment, the time -varying oscillatory x and y components of the 
vortex- core position and their trajectories on the disk plane as 
excited by two (N = 2) Gaussian pulses of field strength H 0 = 
4.6 Oe are plotted in Figs. 3a and 3b, respectively, for AT = 2n/w D 
= 9.09 nsandAT= 1.5 X {2nlco D ) = 13.64 ns. Representative serial 
snapshot images of the XMCD contrast are shown in each second 
panel. In the case of AT = 2tt/cod, where AT is the same as one cycle 
of vortex gyration, orbiting core displacement was increased with 
time upon application of the second pulse. By contrast, for 
AT= 1.5 x (2k/cdd), the core gyration amplitude rapidly decreased 
after application of the second pulse. 

Analytical calculation and numerical simulation. To interpret 
the experimentally observed core gyration amplitude asymmetry 
(amplification versus damping), we conducted micromagnetic 
simulations (see Methods for details). The results are shown in the 
fourth rows of Figs. 3a and 3b, and are in excellent agreement with 
the experimental results, except for slight deviations in their 
oscillation amplitudes and frequencies. These discrepancies can be 
attributed to sample imperfections or several possible simulation/ 
experiment deviations in the dimensions, material parameters and 



field strength. Then, based on Thiele's analytical approach 42 ' 43 , for a 
given magnetization structure taken at a certain time and a given core 
velocity v, all of the forces can be calculated directly from the 
integrations of the magnetization configurations that were 
obtained from the micromagnetic numerical simulations. The 
individual forces are as follows: F G =— G(zxv) with G = 

(L/yM 5 2 ) M-[(dM/dxi) x (dM/dx 2 )]dS, the damping force F D = 



b-Y=^2bjjDj, where a damping tensor — (aL/yM 5 ) 

[(dM/dxj)-(dM/dxj)]dS with the damping constant a, the 



gyromagnetic ratio y, and the saturation magnetization M s . The 
restoring force and Zeeman force are given as Ff +Ff = 

- ^Hf(dMj/dXi)dV with the effective field H eff = SW/SM. 

From these relations, the Zeeman force can be descried as 

¥ H = - ^(SW H /SM)(dM/dXi)dV, where W H = -/x(zxH)-X 

with /i = nRLM s £C(£ = 2/3 for the "side-charge-free" model 35 ). 

Accordingly, the resultant core motion can be determined from 
the numerically estimated individual forces, as represented by the 
magnitudes and directions of all of the forces. For the two cases of 
AT= 1.0 x (2u/wd) and 1.5 x (2tz/wd), temporal variations of the 
force magnitude and orientation are plotted in Figs. 4a and 4b, 
respectively. The definition of the force orientation (0) is illustrated 
in Fig. 4c. The first Gaussian pulse excites the precessional motion of 
the vortex core from its initial position X = 0; if the second Gaussian 
pulse is applied, the Zeeman force F H is generated, and thus, to 
compensate the Zeeman force F H , the magnitude and orientation 
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Figure 3 | (a) Amplification and (b) damping of vortex oscillations. The first row shows the profiles of two serial field pulses of (a) A T = 1 .0 x (2n / cod) 
and (b) AT= 1.5 x (2k/cod). In the second row of (a) and (b), there are serial snapshot XMCD images of the temporal evolution of vortex-core motions. 
The red and yellow symbols in the first row indicate the snapshot time for each image. The third row indicates the x- and /-components and displacement 
of the vortex-core positions (left) from the center position (x, y) = (0, 0) and the constructed vortex-core trajectories (right), as obtained from the time- 
resolved MTXM images; the fourth row shows the corresponding simulation results. 
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Figure 4 | Force-balance relation in vortex oscillation, (a, b) The second 
row presents the strengths of the individual force terms, that is, the gyroforce 
F G , the Zeeman forces F H , the restoring force F w , and the damping forces F D , 
according to the applied Gaussian-pulse fields shown in the first row. The 
bottom row shows the orientations of the individual force terms, that is, the 
angles 6 with respect to the vortex-core-position vector (X), as defined in 
(c). (d) Schematic illustration of orientation and strength of employed 
forces under resonant amplification (left) and forced damping (right), 
where for simplification, damping force F D is not shown, v corresponds to 
the core motion velocity vector. The red dot indicates the core position. 

of gyroforce F G and restoring force F w change. Fig. 4a and 4b show 
the angles and amplitudes of the Zeeman force F H and restoring force 
F w in response to the vortex- core position vector X and applied 
magnetic field H. When the second Gaussian pulse is applied with 
the time interval (AT) of an integer of the vortex precessional period 
1.0 x (2n/w D ), the Zeeman force F H induced by the second pulse is 
in-phase relative to the vortex core orbiting around the center posi- 
tion. To satisfy the force balance, the angle of the gyroforce F G has the 
negative angle, and the magnitude of gyroforce F G and restoring force 
F w increase as shown in Fig. 4a. As a result, the angle of the vortex 
velocity vector v is an acute angle due to the relation of 
F G = — G(z x v). This means that the radius of the vortex gyration 
is increased according to the resultant velocity vector, as shown in the 
first schematic image in Fig. 4d (the damping force F D is not shown 
here for simplicity). In the same manner, the Gaussian-pulse train 
with the time interval of a half integer of the vortex precessional 
period [AT=1.5x (2n/co D )] can reduce the radius of the vortex 
gyration, as shown in the second image in Fig. 4d. 

Furthermore, we experimentally observed the maximum core dis- 
placement |X| max versus AT/AT 0 (AT 0 = 2n/w D = 9.09 ns, 0.5 x 
A T 0 to 2.0 x A T 0 ) by applying a sequence of two Gaussian pulses 
of field strength H 0 = 4.6 Oe and optimal width a = 1 /w D = 1 .45 ns. 
The experimental results (closed circles) of |X| max are compared with 
the analytical calculations (solid line), as shown in Fig. 5. As pre- 
dicted by the analytical form of AT = m x (2n/co D ), the integer 
values of m evidenced amplification of vortex gyration, whereas 
the half integers showed further damping of the vortex gyration. 
This implies that the time at which the second pulse arrives plays a 
role in the in-phase and out-of-phase motions of the vortex core. 

Impact of pulse number on the resonant amplification. Next, we 
experimentally examined the core displacement as a function of time 
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Figure 5 | Maximum vortex-core displacement |X| max for different time 
intervals, AT/AT 0 = 0.5—2.0, where AT 0 = 2n/co D , along with identical 
parameters of pulse number N = 2, o = 1 /co D = 1 .45 ns, and H 0 = 4.6 Oe. 

The yellow dots and solid line correspond to the experimental results and 
analytical calculation, respectively. 



for a different pulse number N (= 1, 2, 3, 4, 5) for the optimal width a 
= 1.45 ns, time interval AT = 9.09 ns, and field strength H 0 = 
4.6 Oe, as plotted in Fig. 6. The inset in Fig. 6 shows the 
corresponding simulation results for the real Py disk and explains 
well how the core displacement increases with the number of 
coherent optimized field pulses. Figure 7a shows the maximum 
displacement of the core position for N = 1, 2, 3, 4, and 5, as 
obtained from Fig. 6, and provides evidence that the amplification 
of the vortex- core gyration by N = 5 is three times greater than that 
by N = 1. In Fig. 7b, for N = 10, the maximum displacement is 
plotted with H 0 , indicating that with application of more Gaussian 
pulses, a sufficient displacement of core motion, 200 nm (1/7.5 of the 
disk radius), can be achieved with the extremely low field strength of 
2 Oe. Interpolation of this data showed that with 1 Oe, 100 nm 
displacement can be achieved. 

Discussion 

We experimentally demonstrated that the resonant excitation and 
amplification of vortex gyrations can be achievable with extremely 
low power consumption as much as ~ 1 Oe by application of a series 
of coherent Gaussian pulses of optimal interval AT = 2n/co D and 
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Figure 6 | Vortex-core displacement as function of time according to 
different number of coherent pulses. The trains of Gaussian magnetic 
field pulses under the optimal conditions o= 1 /cop = 1.45ns, AT = 
1.0 x (2n/co D ) = 9.09 ns, and field strength H 0 = 4.6 Oe were applied. The 
pulse profile according to the time is shown in the first row. The right inset 
shows the corresponding simulation results obtained under the same 
conditions as those prevailing in the experiment. 
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Figure 7 | Experimental results of maximum vortex-core displacement 
versus pulse number in (a) and versus H 0 in (b). The trains of Gaussian 
magnetic field pulses under the optimal conditions a = 1 / Wd = 1 .45ns, A T 
= 1.0 x (2k/cod) = 9.09ns were applied. For comparison, the results in (a) 
also are plotted in (b). 

pulse width a=\/ cdd> Signal gain can be defined as the mean ratio of 
the final vortex-gyration amplitude to the initial vortex gyration. 
Therefore, a large displacement of vortex gyration indicates a suf- 
ficient signal gain detectable through Tunnel magnetoresistance 45 
(TMR) or Giant magnetoresistance 46,47 (GMR) in possible spin value 
structures. Since signal gain in electronic devices is a measure of the 
capacity to increase the power or amplitude of a signal from input to 
output. The present work quantitatively clarifies the fundamental of 
vortex excitations and their amplifications, and might provide a way 
to significantly amplify vortex-gyration-based signal gains in future 
devices at considerably reduced power consumption (~ 1 Oe in field 
strength). 

Methods 

Sample fabrication. For the measurement of soft X-ray transmission through the 
sample, all of the disks and electrodes were prepared on 100-nm-thick silicon nitride 
(Si 3 N 4 ) membranes. The 40 nm-thick Py films were deposited by magnetron 
sputtering under base pressures of less than 5 X 10" 9 Torr. Py disks of radius 2R = 
2.0—6.0 fim were patterned by typical e-beam lithography (Jeol, JBX9300FS) and 
subsequent lift-off processes. Each disk was placed on a strip electrode composed of Ti 
(5 nm)/Cu (70 nm)/Au (5 nm) and patterned, like the Py disks, by e-beam 
lithography and subsequent lift-off processes. 

Time-resolved soft X-ray microscopy measurement. Snapshot images of vortex- 
core oscillations were measured by full-field soft X-ray transmission microscopy at 
the Advanced Light Source (Beamline 6.1.2) in Berkeley, California, USA, using a 
stroboscopic pump-and-probe technique 44 . The microscopy boasts a temporal 
resolution of 70 ps and a lateral resolution of less than 25 nm. The in-plane 
magnetization contrast was provided by X-ray magnetic circular dichroism 
(XMCD) near the Fe L 3 absorption edge (around 707 eV), according to the in- 
plane geometry whereby the sample was mounted at an angle of 60° with respect 
to the propagation direction of incident X-rays. In order to achieve adequate 
XMCD contrasts, 10 to 20 individual images were accumulated. The structural 
contrast was normalized to an image obtained under a saturation field. Based on a 
stroboscopic pump-and-probe technique 44 , Trains of Gaussian- shaped field pulses 
were synchronized to X-ray pulses of ~3 MHz frequency, and the time delay 
between the pulse trains and the incident X-ray was varied from 0 to 60 ns in 
increments of 1.14 ns (eight images per cycle of vortex gyration, at a frequency of 
110 MHz). In the experiment, a Gaussian current pulse was applied with an 
arbitrary waveform generator. The height of the output pulse voltage was 
measured as 350—370 mV with an oscilloscope with an impedance of Z = 
50 Ohm, which corresponds to 0.7~0.74 transmittance. On the basis of Ampere's 
Law, we estimated the resultant Oersted field strength and spatial distribution 
with the measured output voltage and the impedance of the oscilloscope 
employed. Accordingly, the field strength was determined to be H 0 = 4.6 Oe at 
the center of the Py disk (20 nm from the electrode surface). The local distribution 
of the magnetic field was calculated using the analytical equations shown in Ref. 
48. 

Micromagnetic simulations. Micromagnetic simulations were carried out for a 
model system of the same dimensions as those of the sample and under the same 
conditions as those in the experiment. To numerically simulate the vortex-core 
oscillations, we employed the OOMMF code 49 , which utilizes the Landau-Lifshitz- 
Gilbert equation of motion of magnetization:<9M/<9f = — y(M x H e ff) + a/|M| 
(M x dM/dt)) [Ref. 50 and 51], with the gyromagnetic ratio y (2.21 X 10 5 m/As) and 
the phenomenological damping constant a (0.01). We used a unit cell size of 4 X 4 X 



40 nm 3 . To match the eigenfrequency of the model to that of the sample 
{co D = 271 x 1 10 MHz), we used M s = 7.6 X 10 5 A/m for the numerical 
simulations 432 ' 53 . The other simulation parameters for Py were as follows: exchange 
stiffness A ex = 1.3 X 10~ n J/m and magnetocrystalline anisotropy constant K = 0. 
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